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I. INTRODUCTION 



beyond the standard model 



The unification idea, mainly in SU(5) is still an interesting alternative for the physics 



2J. Unfortunately, the minimal non- super symmetric SU(5) 



model has been ruled out by experimental data: i) the proton is more stable than the 
prediction of the minimal model |3j; ii) the value of weak mixing angle at the Z-peak 
sin 2 9w(M z ) = 0.23122(15) [or alternatively a s (M z )] does not agree with experimental 
data [4]. It means that the three coupling constants do not meet at a single point if only the 
standard model particles are taken into account; iii) the electron and d-like quark masses are 
equal at the unification scale, and iv) last but not least, neutrinos are massless in the model. 
Moreover, the supersymmetric version, i.e., SUSY 577(5), although it allows an unification of 
the coupling constants, it has serious problems with the proton decay [5j (however see [(J ) and 
probably also with the electroweak data [7] . Thus, it appears natural to ask ourselves if there 
are other options besides SUSY SU(5) that yield convergence of the couplings, the observed 
value of the weak mixing angle and the other parameters at the Z-pole, an appropriately 
stable proton and, at the same time, realistic fermions masses including neutrino masses. 
Another problem, not necessarily related to the previous one, concerns the existence of 
axions 8J. Recently, the interest in theories involving such particles has raised also due 
experiments devoted to the search of axion-like particles [9]. If the axion does exist it is 
important to know the realistic model in which the Peccei-Quinn (PQ) symmetry can be 
automatically implemented and how the axion parameters can be stabilized against possible 
semi-classical gravitational effects [lol ]. 

On the other hand, it was shown in Ref. ll( that in the context of the multi-Higgs 
extension of the standard model with an invisible axion proposed in Ref. 12] we have: i) 
the unification of the three gauge coupling constants near the PQ scale; ii) the model predicts 
the correct value of the weak mixing angle at the Z-peak; iii) the axion and the nucleon 
are stabilized by the cyclic Z13 <S> Z% discrete symmetries; finally, iv) although neutrinos got 
an arbitrary Dirac mass, through the effective d = 10 operators ApgA -5 L^ u L^ v <p b , the 
left-handed neutrinos get also a Majorana mass < 2 eV and the right-handed neutrinos 
acquire a large Majorana mass term via d = 7 effective operator Apg iy^ R (M R ) ab u bR ((j)*(f)) 2 , 
implementing in this way a see-saw mechanism at the PQ energy scale. 

Here we will consider an 577(5) grand unified theory which unify the model of Ref. 12], 
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in such a way that the partner of the standard model fields in SU(5) multiplets are new 
heavy fields. This model allows an stable proton, unification of the three coupling constants, 
a natural PQ symmetry of the classical Lagrangian, and the axion being protected against 
semi-classical gravity effects. 

The outline of the paper is as follows. In Sec. [TT] we give the representation content 
of the model and the Z13 and PQ charge assignments of the several multiplets. Next, in 
Sec. IHIt we consider the running equations for the three gauge coupling constants related 
to the low energy SU(3)c ® SU{2) L ® U(l)y symmetry group. In Sec [TV] we consider the 
proton stabilization and other phenomenological consequences concerning the model; finally 
the last section is devoted to our conclusions. 



II. NON-SUSY SU(5) GRAND UNIFIED THEORY 



In Ref. [12|] the representation content of the standard model was augmented by adding 
scalar fields and three right-handed neutrinos, in such a way that a discrete Z13 ® Z% sym- 
metry was implemented in the model there. Explicitly, the particle content of the model is 
the following: Ql = (u, d) T L ~ (3, 2, 1/3), Ll = [y, V)\ ~ (1, 2, —1) denote quark and lepton 
doublets, respectively; u R ~ (3,1,4/3), d R ~ (3,1,-2/3), l R ~ (1,1,-2), v R ~ (1,1,0) 
are the right-handed components. It was also assumed that each charged sector gain mass 
from a different scalar doublet: H u , Hd, Hi and H u which generate Dirac masses for u- 
like, c?-like quarks, charged leptons and neutrinos, respectively (all of them of the form 
(1, 2, +1) = (tp + , <£>°) T ). Some other scalar fields were also considered in order to permit the 
full symmetry realization: a neutral complex singlet ~ (1, 1,0), a singly charged singlet 
h + ~ (1, 1, +2) and a triplet T ~ (1, 3, +2). Next, we wonder what is the simplest group 
embedding the above representation content. The answer is: SU(5). To achieve this, along 
with a Z13 symmetry, we have to add new fermions and scalar fields. 

In this vein, the representation content of our SU(5) model is as follows. For each family, 
the fermion representation content under SU(5) D SU(3)c <8> SU{2) L ® U(l)y, there are two 
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where .E and N are heavy charged and neutral leptons, respectively, and U, D are heavy 
quarks having the same electric charge of the respective quarks u, d. Finally, in the fermion 
sector we have to add fermionic neutral singlets {N c )i = iV£ and {v c )l = v c L - We have used 
a notation in which the subindex e(d) denotes the multiplet to which the known leptons (d- 
like quarks) belongs to; on the other hand, the it-like quarks always belongs to the decuplet 
Notice that since the known quarks and leptons belong to different representations of 
577(5), we have to impose that both quarks and leptons, and not quarks and anti- leptons, 
have gauge interactions through the left-handed components 13]. 

The scalars of the model are the usual 24, here denoted by 024, with vacuum expectation 
value (VEV) (024) = M24 diag(l, 1, 1, — §, — |); a complex singlet O which is almost the 
axion (we note that by considering a complex 24 it is possible to implement the axion in 
this model however this may introduce troubles with proton decay). In order to break 
the SU(2) <g> £7(1) symmetry and generate the fermion's Dirac masses we use four Higgs 
multiplets: two 5, and two 45* to avoid the prediction m e (Mjj) = md(Mu) (The using of 



15].) Finally, we add a 10 (-D10) and a 

12]. Ti 5 gives 



45 for avoiding this mass relation was done in Refs. 

15 (T15) which contains, respectively, the singlet h + and the triplet T of Ref. 
Majorana masses to the active neutrinos. We will denote the 5 as H% = {h\, h%, h^, h+, h° a) ) 
with a = e,d\ and their VEVs are (H^ a ) = {v a5 /\/2)8% ; the 45* will be denoted by Hf = 
{H?)f, {Hf)f = {Hf)f = 0, with (Hf)f = (Wv 7 ^ - 4#yfl#; 

finally, (T^f ) = (t> 15 /\/2)#5 5g . The decuplet D w does not necessarily get a VEV at lowest 
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order. Since in this model all scalar's VEVs are of the order of the electroweak scale, except 
024 and 0o which have VEVs of the order of grand unified theory (GUT) and PQ scale, 
respectively, we have still the hierarchy problem. It is only ameliorated because the GUT 
scale is lower (as we will show below) than in other grand unification models. 
Consider the following Yukawa interactions, 



C Y = {^ e ) R [G e ^ eL Hl* + G ei ^ eL Hf + G u v c L Hl] + (^) R eK v ^ eL H^ 



+ (^d) R [G d5 $ dL Hf + G M5 <$> dL Hf + G N N c L H 5 d ] + (H%) R e K d <$> dL Hf* 



- {^ R eF v ^ eL Hl - (^) dR eF d ^ dL H b d - (* e ) R G el5 V c eL T 15 + H.c, (3) 

where G, K, F are 3x3 complex matrices but we have omitted generation and SU(5) indices; 
e denotes the SU(5) fully antisymmetric tensor. With Eq. ([3]) we obtain the mass matrices 
(T denotes the transpose matrix) 

v * v * 
M e = G e5 — 3G e45 w e 45, Mb = G e §-^- + G e ^v e ^, 

M v = V2v e5 (F v + Fy) + 2v / 2< 45 (^ " K v ), (4) 

and 

T ^ d 5 T ^c?5 

Me = G d5 — 3G d45 v d45} M d = G d5 — + G d45 v d 4 5 , 

M u = V2v d5 (F d + Fj) + 2V2v* d45 (Kj - K d ), (5) 

M mrac = ( Ve& /y/2)G%, and M° irac = (v d5 /^2)G T N . The left-handed neutrinos have a 
Majorana mass term coming from the T 15 : M^ Ia ^ orana = (f 15 /v / 2)Gj 15 . Both t> a5 and v a ^ 
are of the order of the electroweak scale, in fact XlaG'^P + l'^sH + l^isI 2 = (246 GeV) 2 , a = 
e,d, with |t>i 5 | < 3.89 GeV [16J. For instance, using only one generation, if t> e5 = f e45 = v e , 
assuming that these VEV are real and neglecting vis, we have, from Eqs. (jlj) and ([51), 
M e = (GJ5/2 - 3G^ 5 )v e and M D = (G e5 /2 + G e45 )v e (and similarly for M E and M d ), so 
we can choose the Yukawa coupling constants must be such that M e <C Mu, M u -C Mu 
and M d <C Me- In the context of three generations all these mass matrices are 3x3 
matrices and those relations among the masses refer to the respective eigenvalues. Right- 
handed components of neutrinos and the neutral leptons N R , get also a Majorana mass term 
through the interactions with the axion \\\ . 

We see that the representation content of the model implies that the vector bosons do 
not induce, at the tree level, the nucleon decay because these interactions involve the usual 



quarks and heavy leptons; or heavy quarks with the usual leptons. The same is true for the 
Yukawa interactions if they are given only by these in Eq. ([3]). This diminish the importance 
of the constraints coming from nucleon decay on the leptoquarks masses. Thus, they may 
have a mass lower than the unification scale. Notwithstanding, when studying the evolution 
of the coupling constants, we will assume that all leptoquarks are heavy enough and do 
not consider them in the running of the couplings. Next, we will show that the Yukawa 
interactions in Eq. ([3]) are the only ones allowed by an appropriate discrete symmetry. 

Let us use the fact that a symmetry with N being a prime number does not have 
any subgroup, in other words, it cannot be decomposed as Z p ® Z q , (p,q < N), so that 
the Z^ symmetry may be a subgroup of a unique local group U(l). In this vain, let us 
introduce the following Z 13 symmetry in the Yukawa interactions in such a way that only 
these interactions in Eq. ([3]) are allowed. The fields of the model transform under Z13 as 
follows: 

{^ C ) eL -> W 3 (^ C ) eL , {^ C ) dL -> 1 (* c ) d£ , $ eL -> Wf^eL, $ dL -> U^dL, 

-^d 5 -> U<sHf, Dio -> ^3^10, ?15 -> 1 Ti5 ! 024 ~> ^0024, 00 ~> ^400- (6) 

We have assumed that the three generations are replicas under Z13. However, it could be 
interesting to consider the case when this is not the case. 

The scalar potential has hermitian quadratic terms /i^x^X (where x denotes any of the 
Higgs scalar multiplets of the model), which are needed to break the electroweak symme- 
try, trilinear and quartic Hermitian terms, and non-Hermitian self-interactions which are 
trilinears 

H e H e 024, H^Hf^i, (Hl) 2 D* m H d H d *-D* , (H d ) 2 D 10 , (7) 

and quartic: 

H d H d Mlfil i^<( 024, H d H d JJ w ( P2A,H d H d H e H e ,H d H d {H d ) , 
(HffT* 5 , 7^02400, {HffDl^ (H 5 e fHt 5 *T l5 , (#f) 2 T 15 0*, {Hi) 2 Hf Hf , 

( Tjb\2 ( Tjib\2 f tt5\2 f tt45\2 tt5 tt4S i rr5 1 2 tj5 tj45\ tt4,5\2 o5 rj-45 I rj-45 1 2 

\ H e) [He ) A H d) [ H d ) ,H e H e \H d \ ,H e H e \H e | ,-H e -H e \M d \ , 

7_x5 7_x45|n 1 2 £j5 £j45 |m 12 rj-5 z_r45 j2 I it5 1 2 r_r5 it45 0-5/ 7^45*^2^ tt5 7_r45 7_x5 rr45 
^e^e l-^iol ,^e^ I J 15 1 ,-n e ^ e 024, l^el ^rf^d , H e[ H e > I 15,H £ H e t± d ti d , 

{HlfHfHf, (tf d 5 m o 024, (^ 5 ) 2 T 1 * 5 0o, HlHf\H?\\ H*Hf\Hf\\ HlHf\D^\\ (8) 
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Moreover, with the interactions in Eq. ([3]) and the non-hermitian interactions in ([7]) and 
(jHJ), allowed by the symmetry in Eq. ([6]), the PQ symmetry is automatic. The PQ charges 
are shown between parenthesis in units of the PQ charge of ^ as follows: 

(M> c ) eL (-l/3), m dL (l), $ £ l(1/9), <M-l/3), i^(5/9), 

iV£(-5/3), HH-2/9), iff (2/9), fl|(2/3), #f (-2/3), 

D 10 (4/3), T 15 (2/3), 24 (O), o (-2/3). (9) 

As in the model of Ref. 12], the Z13 protect the axion against possible semi-classical 
gravity effects. The model has no domain wall problem 

III. EVOLUTION OF THE COUPLING CONSTANTS 

In order to study the running of the coupling constants in a consistent way with the 
present model, we augmented the representation content of the model of Ref. [12j. Hence, we 
assume that the only extra degrees of freedom that are active at low energies, i.e., below the 
unification scale but above the electroweak scale, transforming under the standard model 
(SM) symmetries are (per family): (A, E) T L ~ (3,2,-1), (U, D) T L ~ (3,2,1/3) and the 
respective singlets E R ~ (1,2,-2), N R ~ (1,1,0), U R ~ (3,2,2/3) and D R ~ (3,1,-1/3). 
In the Higgs boson sector we add also four scalar doublets, two of them H^, Hi are those that 
belong to 5 and two others, say H u and H v which belong to the 45; a triplet T belonging to 



T 15 and the singlet h + which is part of the D w . As in Ref. Uj only h + will be considered 
with mass of the order of the unification scale. 

Let us look at the evolution equations at the 1-loop approximation, with all the new 
fermions entering only above an intermediate energy scale fi IS which is certainly bigger than 
the electroweak scale. Some of the new fermions could have mass below the known heavier 
Standard Model particles, but we will not consider such possibility. Below we comment 
more on that (see Sec. |V|). Thus, the 1-loop equations are: 

_L = _J_ _ _L L i n + h ls i n JL\ , (10) 
cn{pi) a t (M z ) 2n[ l M z 1 ^ S J' 1 ; 

where ati(M z ) = gf(M z )/4:7i are the usual gauge couplings defined for these equa- 
tions, and bi the well know coefficients for a general SU(N) gauge group, given by 
bi = (2/3)ET ffi (F) + (1/3) J2T Ri (S) - (11/3) C 2i (G) for Weyl fermions (F) and com- 
plex scalars (5), and T R (I)5 ab = Tr{T a (J), T b (I)} with I = F,S; T R (I) = 1/2 for 



the fundamental representation, C2(G) = N when N > 2, for £7(1), C 2 (V) = 0, and 
TRi(S a , F a ) = (3/5)Tr(y a 2 /4). The same is valid for the b'/ with the counting extending to 
the exotic fermions representations. At the SU(3)c <S> SU(2) L <g> U(l)y energy level with 
N g fermion generations, Nh scalar doublets (Y = ±1) and non-hermitian scalar triplets 
(Y = 2), and N s charged singlets, we have: 

4131 412 22 4 

bi = -N g + —N H + -N T + -N s , b 2 = -N g + -N H + -N T - — , b 3 = -N g - 11, (11) 

where a grand unification normalization factor (3/5) for the hypercharge Y assignment is 
included in b\. So that according the additional representations in the beginning of this 
section only the heavy fermions are activated above fi IS , we have (with N s = 0) 



(h, b 2 , 63) = (5, -2, -7), (b[ S , b 2 , &; S ) = (9, 2, -3) 



(12) 



Note that there is no asymptotic freedom in a 2 at the 1-loop level. The grand unification 
mass scale and the weak mixing angle are given by 



, [fa" 1 (M z ) - fc^" 1 (M z )] 
M GUT = Mrs exp ^ 2tt L5 , T l % r J l_L 1 > x 



(b[ s + f ^ - |C) 



M 7 \ (tf+W-Sbl^ 



"1 T 5 2 5 3 



and 



sin 9 W = - + - 



5a(M z 



8 2tt 



(61 -62) In 



Mz 

Vis 



+ {b 1 -b 



IS , IS \ , / /I, 

in 1 



M, 



GUT 



However, since 



(13) 



(14) 



, , .IS .IS . 6 O IS 6 IS O IS 

b 1 -b 2 = b 1 - b 2 , 61 + -62 - -63 = 6 X + -o 2 - -63 , 

5 5 5 5 



(15) 



Mqut, and sin 2 Ow(Mz) defined at the Z boson mass, do not depend on the scale fi IS and 
we are left with 

a' 1 (M z ) - l^ 1 (M z y 



M GUT = M z exp 



2tt- 



|6i + b 2 - §6 3 



and 



sin 2 



5a(M z 



8 2tt 



(61-62)111 



GUT 



(16) 



(17) 



Both Mqut and sin 2 ^vy(M^), at the one loop level, are the same as in Ref. Us- 
ing M z = 91.1876 GeV; a(M z ) = 1/128; and a 3 (M z ) = 0.1176 Q, we obtain 
Mqjjx = 2.86 x 10 13 GeV and sin 2 ^(M z ) = 0.23100, in agreement with the usual 
value [4]. Moreover, using the evolution equations in Eq. ffTOl we get a GUT — 23(21), if 
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FIG. 1: In this figure we show the convergence point for the SU(5) model here. 



fj, IS ~ 1 TeV (fi IS = Mz). The inclusion of the scalar singlet at low energies (N s = 1) gives 
worse values for this mixing angle so, it must be considered with mass near the unification 
scale. In Fig. 1 we show the evolution of the coupling constants at the 1-loop level in the 
present model. An analysis at the 2-loop can be done, but in general it does not lead to a 
prediction unless the top quark and all extra fermion and scalar fields are taken into account. 
This result in a large set of coupled equations [3] that deserve more careful study. 

Notice that in this extension of the SM, the Yukawa interactions can be similar to those 
in Ref. 12] but, it is worth noting that, if we want to avoid a general mixing in each charge 



sector the extra quark (lepton) generation must transform under Z13 in a different way from 
those of the usual lepton (quarks). However, we recall that getting a small mixing it can be 
interesting if in the future a departure from unitarity in the Cabibbo-Kobayashi-Maskawa 
mixing matrix would be observed [^(J. In fact, even the usual three generations can be 
transformed under Z13 different from each other in such a way that predictive mass matrices 
can be obtained. 



IV. STABILIZING THE PROTON 



As we said before, in this model nucleon decays are forbidden at the tree level. Here 
we will discuss this point in more detail. The effective operators with dimension six 21], 
d = 6, that can induce the proton disintegration do not operate in our model because vector 



leptoquarks always mix the usual fermions with the heavy ones. However, without the Z%s 
symmetry, there are still dangerous d=4 operators coming from the Yukawa couplings with 
the H A5 Higgs scalars. For instance, without that discrete symmetry, Yukawa interactions 
like ^en^dLH^ 5 and $^$^eifg 5 * are allowed. These terms induce the proton decay through 
interactions like Q c mR ^o ■ fj m L L , and e mnp Q c mTi ea ■ ff n Q pL , respectively; here e = ia 2 , and 
m,n,p, are color indices, and ff m is the colored scalar triplet belonging to H^ 5 . Once the 
Zi 3 symmetry is introduced the Yukawa interactions allowed are just those given in Eq. 
and they only induce interactions like Q c mR ^ & • rj m L' L and emn P QmR e & ' VnQpLi where the 
primed fields are heavy quarks, U, D, or heavy leptons, E, N. Hence, with the interactions 
in Eq. ([3]), independently of the mixing in the scalar sector, the nucleon is not allowed to 
decay at the tree level. The model is in this respect phenomenologically safe. 

The Z13 symmetry introduced in this model allows effective interactions with flavor chang- 
ing neutral current. For instance, 



here a, b, c, d are family indices, and gz N and M Zn denote the coupling constant of the 
Z N C U(l) i oca i and the mass of the (heavy) vector boson associated with this symmetry, re- 
spectively and habcd are dimensionless constants (there are also effective interactions induced 
by the heavy scalar that condensate at very high energies). The interactions in Eq. (118p 
induce rare transitions like \x — > eee. Neglecting the electron masses we can write the width 
of this decay in terms of the muon decay width as follows: 



where g 2 and Mw are the well-known parameters of the standard model and we see that 
even if gz N ~ 0(g) with Mz N > 10 3 Mw, we have already got a suppression factor of 1CT 12 . 
However, it is more natural that Mz N be of the order of the breakdown of the local U(l) 
symmetry i.e., at least of the order of the PQ scale. It may be also interesting to assume 
that gz N «C g at low energies, in such a way that, since gz N which is not an asymptotic 
free parameter it can fit with g at a high energy and the other coupling constant of the low 
energy model. 

Just as another example, there are also interactions like 




^— h abcd L aL ^L bL L cL 7 M L di , 



(18) 




(19) 




^KbcdQaLYQ b LF cL ^F dL , 



(20) 
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where F = Q,L and hi is another dimensionless matrix. When F = Q this interaction will 
induce a contribution to AMk and other related parameters. Notice, however, that 

AM K oc (°*L G F B K f K m K , (21) 

V 9 M z N J 

we see that this contribution to AMk is rather small for the same values of the param- 
eters in Eq. ( |T9l) . 

In general discrete symmetries may be not free of anomalies. Although it is interesting to 
looking for cyclic local discrete symmetries that are anomaly free, we would like to emphasize 
that it is not necessarily a loophole of models with anomalous Zn symmetries. If gz N < g 
the transition violating B and L conservation induced by the anomaly of the Z^ symmetry 



will be smaller than e ~ 16 ^ 2 /9 2 « lQ- n7 S 



221 ]. with £ = 0(1) a model dependent parameter. 



Although this transition is negligible at zero temperature it may be important in earlier ages 
of the universe as a mechanism for leptogenesis generation through the decays of the heavy 
neutral leptons if CP violation is implemented in it. In fact, the model allows several ways 



to implement baryogenesis and leptogenesis 23|, |24( as we will shown elsewhere. 



V. CONCLUSIONS 



Summarizing, we have obtained an SU(5) extension of a previous model of Refs. 12], 
which is as good as SUSY SU{5) with respect to the unification of the electroweak and strong 
interactions. We also have in this model that the proton is stable and the PQ an automatic 
symmetry of classical Lagrangian, with the invisible axion protected against gravitational 
effects by a local Z 13 symmetry. The unification of the three coupling constants occurs at the 
PQ scale as in 11|. It is important to realize that the model does not admit supersymmetry 
at least at low energies, but we have seen that, in order to have unification, supersymmetry 
is not an indispensable factor anymore. As it was put forward in Ref. the PQ energy 
scale can be related with the mass of the sterile neutrinos and in the present model 
the PQ energy scale is related with the GUT scale. Put all this together and we have that 
in the present model it is possible to have M VR -C M P q ~ M GUT or M UR ~ M P q ~ M GUT , 
depending on the Zis charges assignment. 

Although the scalar and vector leptoquarks do not induce the nucleon decay at the tree 
level, their masses are of the order of the unification scale (they gain masses from the (24)). 
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Since this scale is of the order of the PQ scale it means that both energy scale can be related 
to each other. On the other hand, the exotic quarks and leptons U, D, E gain mass from 5 
and 45 which have VEVs of the order of the electroweak energy scale and, for this reason, 
they could not be very heavy. We recall that the experimental limits on the exotic leptons 
and quarks like E, N and U, D, respectively, are model dependent but since all of them gain 
mass from VEVs of the order of the electroweak scale they must not be very heavy indeed. 
For instance, from data we have lower bounds on the masses (in GeV) of a possible fourth 
family 4|: for sequential E^ 1 charged lepton, we have m > 100.8, CL=95% (decay to uW); 
for stable charged heavy leptons m > 102.6, C.L.=95%; for stable neutral heavy lepton the 
limits are m > 45.0, C.L.=95% (Dirac) and m > 39.5, C.L.=95% (Majorana). Finally, 
for extra quarks of the 6-type (b' 4th generation) the lower limits are m > 190, C.L.=95% 
(quasi-stable b') or m > 199, C.L.=95% (neutral currents); if it decays in 11 +jets, Z+jets, we 
have m > 128, C.L.=95%. Of course, these limits are strongly model dependent (in some 



models the fourth family is almost degenerate 25J]). On the other hand, as we said before, 



irus 



the model has a general mixing among the fields of the same electric charge sector, t 
the generalized unitarity triangle analysis of the Cabibbo-Kobayashi-Maskawa matrix 2(3] 
can be used for deriving upper bounds on the coefficients of the effective operators inducing 



such mixings 261 ] . 



The model has right-handed neutrinos, so it is possible that an 5O(10) would be more 

n 

appropriate for the unification of the model of Ref. [121 ] . However since SU(o) C 5*0(10) 
our SU(5) model is already good enough for implementing a GUT theory for the extension 
of the standard model with proton stability. 
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